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Abstract 

We present results from a jet energy loss calculation using the Gyulassy-Levai-Vitev (GLV) formalism and bulk medium evolution 
from the covariant transport MPC. At both RHIC and LHC energies we find that realistic transverse expansion strongly reduces 
elliptic flow at high p T compared to calculations with transversely "frozen" profiles. We argue that this is a generic feature of GLV 
energy loss. Transverse expansion also leads to stronger high-/? T suppression, while fluctuations in energy loss with the location 
of scattering centers weaken the suppression. But, unlike the reduction of v 2 , these effects nearly disappear once a s is adjusted to 
reproduce R AA in central collisions. 
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1. Introduction 

Understanding parton energy loss in ultrarelativistic heavy-ion reactions has been the focus of considerable recent 
theoretical effort. A variety of phenomenological approaches (e.g., |Q] |2] [3|) formulate the problem in terms of a 
local energy loss rate dE/dL = - f(E(L), T(L), L) along the Eikonal parton trajectory, given by the local temperature, 
position, and parton energy. In the small-coupling regime, more rigorous treatment is possible based on perturbative 
QCD[4, 5, 6|. This includes quantum interference effects and also fluctuations, namely, energy loss along a given jet 
trajectory becomes a stochastic variable that is in general a function of the scattering and emission history of the jet. 

A critical step in computing heavy-ion observables from any energy loss model is spatial and temporal averaging 
over the bulk medium formed in the collision. We employ here the Gyulassy-Levai-Vitev (GLV) framework|6| in 
which a high-energy parton loses energy through gluon radiation induced by interactions with static Yukawa scatterers 
in the medium. It is natural to combine this approach with parton transport for the bulk evolution, such as Molnar's 
Parton Cascade[7| (MPC), because it directly provides scattering center information. 

Our approach is similar to recent work by Buzzatti and Gyulassy 1 8 1, but with a few key differences. Unlike 1 8 1, at 
present we only focus on light partons, and do not include multiple gluon radiation, elastic energy loss, or energy loss 
fluctuations due to variations in radiated gluon momentum. However, we do include realistic 3D medium evolution 
with both longitudinal and transverse expansion, which turns out to influence energy loss and, especially, elliptic flow. 



2. Radiative energy loss and medium averaging 



We consider here the leading n — 1 (single scattering) term in the GLV opacity expansion of the radiated gluon 
spectrum! 6 1 
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Figure 1. Neutral pion suppression factor Raa at midrapidity in central and mid-central Au+Au at tJsnn = 200 GeV from calculations with GLV 
energy loss. Left: results for a medium undergoing longitudinal Bjorken expansion only. Right: results for realistic medium evolution with both 
longitudinal and transverse expansion (see text for details). In both cases, two scenarios are considered, one based on the average energy loss along 
the jet pass (solid), and one utilizing stochastic energy loss AE(z) (dashed). Data from PHENIX (boxes) are also shown| 12 | to guide the eye. 



where the original hard scattering is at z = 0, p(z) is the local Debye screening mass, cr = 9na 2 / (2fj. 2 ) is the (screened) 
total gg — > gg scattering cross section, and x ~ J dzp,cr is the opacity. We integrated this spectrum numerically 

with kinematic bounds k < xE, q < V6£T, and xE > fi to obtain a momentum-averaged energy loss AE m (z) — 
J dxd 2 kEx(dN w /dxd 2 k) for fixed z, i.e., retained energy loss fluctuations due to variations in z only. The probability 
for the scattering to occur at z is p(z) - p(z)cr(z)/x, so the fully averaged energy loss is AE m — J dz p(z) A£ (1) (z). 
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Figure 2. Same as Fig.JTjbut with a slight tuning of a s in each of the four scenarios to match the data above pj > 8 GeV in central collisions. After 
the tuning, all scenarios give practically the same result. 



As customary, in non-static media we reinterpret p(z) in the GLV formula as p(z, t — to + z) along the parton 
trajectory. The density evolution was obtained from the parton transport MPC, employing 2 — > 2 interactions for 
massless gluons. The scattering rate was adjusted to generate substantial V2(pr ~ 3 GeV) ~ 0.25 in collisions 
with b = 8 fm impact parameter, and we set growing o~ gg ^ gg oc t 2 ^ to keep the shear viscosity to entropy ratio 
approximately constantj9| . Initial conditions for Au+Au at ^JJnn = 200 GeV and Pb+Pb at y/s NN = 2.76 TeV were 
based on diffuse Woods-Saxon nuclei. For the transverse density, binary collision profiles were used, while the impact 
parameter dependence of rapidity densities dN{b)/dy was proportional to N par , with dN(Q)/dy = 1100 (Au+Au) 
and 2400 (Pb+Pb) to match the observed charged particle yields. Because we are only interested in observables at 
midrapidity, we set up boost invariant conditions in the coordinate rapidity window \rj\ < 5, with formation time 
t = 0.6 fm. 
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Figure 3. Differential no elliptic flow V2(pr) at midrapidity in central and mid-central Au+Au at -Jltf/i = 200 GeV from calculations with GLV 
energy loss. The value of or 4 was tuned to reproduce Raa in central collisions (cf. Fig. [5). Left: results for a medium undergoing longitudinal 
Bjorken expansion only. Right: results for realistic medium evolution with both longitidunal and transverse expansion (see text for details). In both 
cases, two scenarios are considered, one is based on the average energy loss along the jet pass (solid), and one utilizes stochastic energy loss AE(z) 
(dashed). With realistic 3D medium evolution we find a striking 40-50% reduction in elliptic flow compared to the case of simplified ID dynamics. 
Data from PHENIX (boxes) are also shown| 13 1 to guide the eye. 



Using tabulated densities from the transport, we set the local temperature assuming a massless gas of gluons 
p x 2T 3 and the Debye mass via ft = gT w 2T. At early times r < to we assume linear density build-up[8| 
p = tp(xt,to)/to. We roughly account for additional energy loss off dynamical (recoiling) scattering centers! 10 1 
[(q 2 +fi 2 ) 2 — > q 2 (q 2 +fi 2 ) in ([I}] and elastic scattering lfTTI viarescaling our opacities ;f — > ^/Z with Z = 0.35. Initial jet 
momentum distributions in p+p, Au+Au and Pb+Pb were computed from leading-order (LO) perturbative QCD with 
one-loop running coupling a s (Q 2 ), using CTEQ5L parton distribution function parameterizations with Q 2 = p\ panon - 
Nuclear effects such as shadowing were ignored but isospin (proton-neutron difference) was included. After energy 
loss, jets were fragmented independently using LO BKK95 fragmentation function parameterizations with scale factor 
Q 2 = p^, j, and we assumed 7Tq = {n + + n~)/2 for the neutral pion yield. This procedure reproduces high-pr ttq and 
charged particle spectra in p+p at RHIC and LHC with modest K^lo ~ 2.5 to account for higher-order contributions. 

Below we focus on two basic high-pr observables for neutral pions at midrapidity, the nuclear suppression factor 
Raa and the momentum anisotropy (elliptic flow) Vi = (cos 2</>) PT . Only energy loss was considered, i.e., contributions 
by the radiated gluons to the final spectrum and feedback on the bulk medium due to the jet were ignored. 



3. Main results 

We considered four scenarios based on i) whether the medium is only undergoing Bjorken expansion ("ID" as 
in |8|) or transverse expansion as well ("3D"); and ii) whether average energy loss is used or the stochastic AE(z). 
Figure [T] shows our results for Raa at RHIC, for the same a s = 0.29. In the stochastic case energy loss effects are 
noticeably weaker, which is natural for convex parton spectra ("curving up" at high p T ). We also find that realistic 
transverse expansion significantly enhances jet quenching, which is a generic GLV feature coming from the interfer- 
ence term in ([TJ. Scatterings at large z induce larger energy loss, and with a transversely expanding density profile 
there is higher chance to scatter further away from the production point than in the transversely static case. 

Unfortunately, without precise control over a s , Raa alone cannot differentiate between these four scenarios. As 
shown in Fig. [2] after a slight tuning of a s to reproduce the suppression in central collisions, differences in Raa largely 
disappear. On the other hand, striking difference in v>2 remains between ID and 3D evolution at both RHIC and LHC 
energies, as shown in Figs.[3]and|4] The strong 40-50% reduction of vo at high-p^ with realistic transverse expansion 
is another generic consequence of interference in GLV. Scattering points that lead to most energy loss are biased to 
occur away from the production point and so later in time, by when the expansion makes the system more cylindrical, 
reducing the spatial azimuthal asymmetry that drives elliptic flow. We expect that this strong effect will be manifest 
in more full-fledged GLV calculations as well, such as [8 1. 
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Figure 4. Same as Fig. [3] but for Pb+Pb at 

■\Js~nn — 2.76 TeV and with a s tuned to reproduce Raa « 0. 18 at pj = 6 GeV in central collisions. 
With realistic 3D medium evolution (right) we only calculated \>i for non-central collisions. Transverse expansion reduces V2 dramatically at LHC 
energies as well. Charged particle \'i(pt) data from CMS (boxes) are also shown| 14] to guide the eye. 



At the conference we also presented results based on scattering center ensembles from the transport (not just 
density information) but due to page limitations these will be written up elsewhere. 

4. Conclusions 

We investigated GLV energy loss using bulk medium evolution data from the covariant transport MPC. We find that 
realistic transverse expansion strongly suppresses elliptic flow at high p T compared to calculations with transversely 
"frozen" profiles (as in JSJ). We argue that this is a generic feature of GLV energy loss, raising the question whether 
GLV produces too little elliptic flow at high pj. 

Transverse expansion also enhances the high-p^ suppression, while fluctuations in energy loss with the location of 
scattering centers reduce energy loss effects. However, unlike for V2, these effects nearly disappear once calculations 
are adjusted to reproduce Raa in central collisions. 
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